ارائه مدل عددی گالرکین ناپیوسته IMPES برای مدلسازی آلاینده های زیر زمینی امتزاج ناپذیر با کمک روش Lax-Wendroff
Authors
Abstract:
سابقه و هدف مدلسازی عددی جریان های امتزاج ناپذیر در محیط متخلخل از جمله مباحثی است که بدلیل کاربرد آنها در پایش انتقال آلاینده ها، حرکت آب و نفت در مخازن نفت و علوم هیدرولوژی همواره مورد توجه محققین قرار می گیرد. در این تحقیق، به ارائه یک مدل عددی دوبعدی گالرکین ناپیوسته جریانهای امتزاج ناپذیر در محیط متخلخل با استفاده از استراتژی حل معادلات فشار-ضمنی درجه اشباع صریح (IMPES) مرتبه بالا پرداخته شده است. در معادلات مورد استفاده متغیر های اصلی فشار و درجه اشباع آب می باشند. در این ترکیب عددی برای اولین بار از روش لاکس- وندروف مرتبه دوم در حل معادله درجه اشباع آب استفاده شده است که بعنوان نوآوری اصلی این مقاله تلقی می گردد. مواد وروش ها به منظور مدلسازی عددی آلاینده های زیرزمینی امتزاج ناپذیر، از گسسته سازی مکانی دارای بقای محلی گالرکین ناپیوسته استفاده شده است. برای گسسته سازی زمانی معادله بقای جرم و درجه اشباع (انتقال) آلاینده نیز به ترتیب از روش های اولر ضمنی مرتبه اول و روش مرتبه بالای لاکس-وندروف صریح مرتبه دوم بهره برده شده است.همچنین به منظور بهبود نتایج در تسخیر شوکها و محل ناهمگنیها از تثبیت شارهای تبادلی و نگاشت میدان سرعت در فضای برداری H(div) استفاده شده است. در انتهای هر گام زمانی نیز نوسانات مقادیر درجه اشباع با استفاده از محدودکننده شیب چاونت-جافر اصلاح شده حذف و نتایج تثبیت می گردند. یافته ها روش مرتبه دوم لکس-وندروف بر مبنای بسط تیلور و ترمهای مرتبه بالای مشتق زمانی، دارای دقت قابل رقابت با روشهای مرسوم در استراتژی IMPES همچون روش چند مرحله ای رانج-کوتا گالرکین ناپیوسته (RKDG) بوده و هزینه محاسبات کمتری نسبت به روشهای چند گامی دارد. هر چند اندازه گامهای زمانی و عدد کورانت با توجه به حل صریح معادله درجه اشباع در این روش دارای محدودیت هایی می باشد. نتیجه گیری صحت سنجی مدل تهیه شده با استفاده از مسئله بنجمارک لورت باکلی ارزیابی شده و نتایج حاصل از مدلسازی با نتایج سایر محققین مقایسه گردیده و تطابق مطلوبی بین آنها حاصل شده است. همچنین ارزیابی کارایی و توانمندی مدل با کمک مسائل نمونه برای آبخوانهای بسیار ناهمگن بررسی شده است. نتایج بیانگر آنست که بعلت استفاده از روش گالرکین ناپیوسته دارای بقای محلی و تکنیک های تثبیت کننده شارهای تبادلی، وضوح نتایج مطلوب بوده و محل تماس دوفاز امتزاج ناپذیر بخوبی تسخیر شده و پخش عددی مشاهده نمی گردد. در این مدلسازی از مقادیر پنالتی 50 تا 100 برای نسخه SWIP استفاده شده است که با توجه با مقیاس نمودن ترمهای پنالتی مدل حساسیت چندانی نسبت به بزرگی آن ندارد ولیکن در نسخه OBB مقدار پنالتی صفر می باشد. در انتها نیز آنالیز حساسیت مدل به ازای تغییرات پارامترهای موثر در مدلسازی مدنظر قرار گرفته شده است.
similar resources
توسعه یک مدل عددی جدید بر مبنای گالرکین ناپیوسته برای شبیه سازی تهاجم آب شور دریا به آبخوان های ساحلی
سابقه و هدف: آبخوان های ساحلی از مهم ترین منابع تامین آب شیرین در بسیاری از کشور های جهان، بخصوص در نواحی خشک و نیمه خشک به شمار می رود. بدلیل مجاورت و ارتباط آبخوان های ساحلی با آب شور دریا و تهدید ناشی از آلوده شدن آن ها بواسطه پیش روی آب شور، مدیریت و حفاظت این منابع آب شیرین ساحلی امری کاملاً ضروری است. بنابراین هدف از پژوهش حاضر، توسعه یک مدل عددی جدید برای شبیه سای انتقال آلودگی در آبخوانه...
full textحل عددی معادلات دیفرانسیل معمولی کسری با روش گالرکین ناپیوسته موضعی
در این مقاله، روش گالرکین ناپیوستهی موضعی برای حل معادلات دیفرانسیل معمولی با مرتبهی کسری را در حالت کلی به کار میبریم. در این روش انتخاب (طبیعی) شار عددی آپویند، ما را قادر میسازد تا مسائل مقدار اولیه برای معادلات کسری معمولی را به صورت بازه به بازه و پیشرو در زمان حل کنیم. این بدین معنی است که ما بایستی در هر زیربازه به حل یک دستگاه معادلات از مرتبه پایین $(k+1)times (k+1)$...
full textمطالعه عددی اندرکنش تونل های زیر زمینی با گسل نرمال
در بسیاری از موارد، تونل ها به عنوان زیرساخت های حیاتی در حمل و نقل شهری در مناطق لرزه خیز و در مجاورت گسل ها احداث می شوند. در هنگام وقوع زلزله، در اثر جابجایی تفاضلی طرفین گسل، سطح زمین دچار گسیختگی شده و می تواند به سازه های زیرسطحی در مجاورت گسل خسارات جدی وارد سازد. بررسی زلزله های گذشته مانند زلزله های ترکیه و تایوان در سال 1999 نشان می دهد که با اجتناب از ساخت و ساز در پهنه گسل نمی توان ...
full textPositivity-Preserving Discontinuous Galerkin Methods with Lax-Wendroff Time Discretizations
This work introduces a single-stage, single-step method for the compressible Euler equations that is provably positivitypreserving and can be applied on both Cartesian and unstructured meshes. This method is the first case of a singlestage, single-step method that is simultaneously high-order, positivity-preserving, and operates on unstructured meshes. Time-stepping is accomplished via the Lax-...
full textThe Near - Stability of the Lax - Wendroff Method
In discussing finite difference methods for the solution of hyperbolic par t ia l differential equations, STETrER [1] used est imates on some absolutely convergent Fourier series to prove stabi l i ty and instabi l i ty with respect to uniform convergence. I f / , a complex valued function on the circle, has an absolutely convergent Fourier series, then the n-th power of ] also has an absolutel...
full textMy Resources
Journal title
volume 26 issue 2
pages 1- 27
publication date 2019-06-22
By following a journal you will be notified via email when a new issue of this journal is published.
Hosted on Doprax cloud platform doprax.com
copyright © 2015-2023